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Abstract 

Aims. 

We study the forced rotation of Titan seen as a rigid body at the equilibrium Cassini state, involving the spin-orbit 

sy nchronizat ion . 

Methods. 

We used both the analytical and the numerical ways. We analytically determined the equilibrium positions and the 
frequencies of the 3 free librations around it, while a numerical integration associated to frequency analysis gave us a 
more synthetic, complete theory, where the free solution split from the forced one. 
Results. 

We find a mean obliquity of 2.2 arcmin and the fundamental frequencies of the free librations of about 2.0977, 167.4883, 
and 306.3360 years. Moreover, we bring out the main role played by Titan's inclination on its rotation, and we suspect 
a likely resonance involving Titan's wobble. 
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1. Introduction 

■ Since the terrestrial observations of Lemmon et al. (jl993p . 
, the rotation of Titan, Saturn's main satellite, has been 

■ assumed to be synchronous or nearly synchronous. This 
' has been confirmed by Lemmon et al. (|1995ap and by 
, Richardson et al. (|2004p with the help of Voyager I im- 
' ages. In this last work. Titan's rotation period is esti- 
' mated at 15.9458 ± 0.0016 days, whereas its orbital period 
' is 15.945421 ± 0.000005 days. 

[ The spin-orbit synchronization of a natural satellite is 
very common in the solar system (such as for the Moon and 
the Galilean satellites of Jupiter) and is known as a Cassini 
state. This is an equilibrium state that has probably been 
reached after a deceleration of the spin of the involved body 
under dissipative effects, like tides. 

Recently, Henrard and Schwanen (2004') have given a 
3-dimensional elaborated analytical model of the forced ro- 
tation of synchronous triaxial bodies, after studying the 
librations around the Cassini state. This model has been 
successfully applied by Henrar d on the Galilean satellites 
lo (2005a) and Europa (2005b), seen as rigid bodies. Such 
studies require knowing some parameters of the gravita- 
tional field of the involved bodies, which cannot be con- 
sidered as spheres. Another analytical study has been per- 
formed for Mercury by D'Hoedt and Lemaitre l|2004p . for 
the case of a 3 : 2 spin-orbit resonance. 
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Since the first fly- bys of Titan by the Cassini spacecraft, 
we have a first estimation of the useful parameters, more 
particularly Titan's J2 and C22 (Tortora et al. I2006p . so 
a similar study of Titan's rotation can be made. In this 
paper, we propose a study of Titan's forced rotation, where 
Titan is seen as a rigid body. The originality of this study 
over Henrard's previous studies is that we use both the 
analytical and the numerical tools and compare our results. 

2. Expressing the problem 

Titan is here considered as a triaxial rigid body whose prin- 
cipal moments of inertia are written respectively as A, B, 
and C, with A<B <C. 

2.1. The variables 

Our variables and equations have already been used in 
previous studies; see for instance Henrard and Schwanen 
(|2004p for the general case of synchronous satellites, 
Henrard (|2005ap for lo, and Henrard (|2005bp for Europa. 

We consider 3 reference frames : the first (ei, 62, 63) is 
centered on Titan's mass barycenter and is in translation 
with the inertial reference frame used to describe the orbital 
motion of the Saturnian satellites in the TASS1.6 theory 
(see Vienne & Duriez I1995P . This is a cartesian coordinate 
system whose origin is the center of Saturn, and it refers to 
the equatorial plane of Saturn and the node of this plane 
with the ecliptic at J2000. The second frame (ni, ^2, 't-s) 
is linked to Titan's angular momentum, and the third one 
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if It f 2, fa) is rigidly linked to Titan. In this last frame, 
Titan's matrix of inertia is written as 



/ A ■ 
1= i B 
\0 C 



(1) 



We first use Andoyer's variables (see Andover 119261 and 
Deprit fT967p . which are based on two linked sets of Euler's 
angles. The first set {h,K,g) locates the position of the 
angular momentum in the first frame (ei, 62, 63), while the 
second {g,J,l) locates the body frame (/i,/2,/3) in the 
second frame tied to the angular momentum (see Fig.[T]). 

The canonical set of Andoyer's variables consists of the 
three angular variables l,g,h and their conjugated mo- 
menta defined by the norm G of the angular momentum 
and two of its projections: 



/ 

9 

h 



L ^ G cos J 
G 

H = GcosK 



Unfortunately, these variables present two singularities: 
when J = (i.e., the angular momentum is colinear to /a, 
there is no wobble), I and g are undefined, and when K = Q 
(i.e., when Titan's principal axis of inertia is perpendicular 
to its orbital plane), h and g are undefined. That is why we 
use the modified Andoyer's variables: 



p = I + g + h 
r = —h 

^9 = 




P{l-coaK) 



where n is Titan's mean orbital motion , q = —I, and Q = 
G—L = G(l — cos J) = 2G'sin^ ^. With these new variables, 
the singularity on I has been dropped. 



2.2. The free rotation 

To describe the dynamics of the system, we should con- 
sider the free rotation and the perturbations by other bod- 
ies. The Hamiltonian of the free body rotation is also the 
kinetic energy of the rotation T — i(a;|G) where a; is the 
instantaneous rotation vector, G the angular momentum 
vector with respect to the center of mass, and {(j^\G) the 
scalar product of the vector u) and G, where uj and G are 
respectively defined as 



(2) 



and 



and consequently 



sin^ / 



B 



7^2 



As a result, the Hamiltonian of the free rotation in the 
modified Andoyer's variables is 



T 



nP 
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71+72 ^2 71-72 2 



1-71-72 



7i +72 



with 



and 



2C -A- B MR^ 
71 = — = J2- 



2C 



C 



2C 



2.3. Perturbation by Saturn 



B-A MR' 
72 = = ^022- 



C 



(7) 
(8) 

(9) 



Considering the parent body Saturn as a point mass M-^, 

the gravitational potential of the perturbation can be writ- 
ten as 



V 



pdW 



(10) 



w 



where p is the density inside the volume W of the body 
and d! the distance between the point mass and a volume 
element inside the body. Using the usual expansion of the 
potential in spherical harmonics (see for instance Bertotti 
and Farinella ll990p . we find 



V = 



n 



[C™ cos milj + 5™ sin mtp] 



(11) 



where ^ and (f) are respectively the longitude and the lati- 
tude of Saturn's barycenter of mass in Titan's frame, and d 
the distance between this Saturn's barycenter of mass and 
the origin of the frame (Titan's barycenter of mass). If we 
limit the expansion of ()11|) to the second order terms and 



drop the term — j^, which does not produce any effect on 
the rotation, we have 



G = Aojih+Boj2h + CuJih. 



(3) 



We also deduce from the definitions of the angles I and 
J (the wobble): 



V = -- 



2^3 



^MR' 



J2{x^ +y^) + 2c22ix' -y^) (12) 



G = G sin J sin ^/i + G sin J cos Z/2 
from which we can easily deduce 



G 
A 



sin J sin Ifi 



G 
B 



where x, y, and z are the coordinates of Saturn's cen- 
ter of mass in Titan's frame (/i, /2, /s) (so we have d = 
G cos J fs , (4) a/x^ -I- y2 _|_ 2;2) . Here, d depends on the time since Titan's 
motion around Saturn is not circular, but we can introduce 
do, the mean value of d, since oq is Saturn's mean semi- 
major axis and cq its mean eccentricity. (They correspond 
Q respectively to Titan's mean semimajor axis and mean ec- 

centricity in a Saturnian frame.) 



sin J cos ;/2 + — cos J/3 (5) 
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We use the formula 



with 



do = a{ 1 + — 



(13) 



coming from the development of - (see Brouwer & 
Clemence[19611): 



where R-^ is Saturn's radius, and J^f^ its J2. 

Finally, the Hamiltonian of the problem reads 



(20) 



- = 1 + — + -e+-e^ cosTW 
a 2 \ " / 

p2 Q 

-— COS27W - -e^cosSM + 0(e^) 
2 8 



(14) 



from which the mean anomaly A4 disappears after averag- 
ing. The perturbing potential V now reads 



and we set n*^ — ' . 

"0 



so that we can write 



V 



do 



with 



and 



Si 



3 /n* 



2 V 71 



71 



(15) 



(16) 



(17) 



(18) 



where M and R are respectively Titan's mass and radius. 

As Henrard (j2005bp did for Jupiter, we also take 
Saturn's oblateness into account. The perturbing potential 
due to Saturn's oblateness reads 



V, = 6sCn^ (^^) ' [Si{x' + 2/2) + S2ix' - y')] (19) 



+n 



nP 



do 



71 + 72 p2 



n 

71 - 72 

1 - 71 + 72 




1 - 71 - 72 

' [5i{x^ + y^) + 52{x'-y^)\. 



(21) 



3. Analytical study 



We intend to use the Hamiltonian (|2T|) to analytically deter- 
mine the equilibrium position of Titan in the Cassini state 
related to the spin-orbit synchronization and the 3 frequen- 
cies of the free librations around this equilibrium, using 
the method explained in (Henrard & Schwanen I2004p . For 
this analytical study, we consider that Titan has a circular 
orbit around Saturn, whose inclination on Saturn's equa- 
torial plane is given by only one periodic term extracted 
from TASS1.6 ephemerides (Vienne & Duriez |1995i This 



implies that the ascending node of Titan oscillates around 
a fixed value, so it cannot disappear after averaging the 
equations. That is why the analytical solutions of (Henrard 
and Schwanen I2004p cannot be used directly, and we first 
must check that they become unchanged without averaging 
the ascending node. The true orbital eccentricity of Titan 
is about 0.0289, but the opportunity to neglect it will be 
discussed later, after comparison with the numerical study. 

This way, the vector locating Saturn's barycenter is co- 
linear to 



XiB-i +yie2 + ZiEs 



(22) 
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with 



- ( cos ^6 cos(A6 - fZe) - cos h sin sin(A6 - fZe 



(23) 



2/i = " ( sin Jig cos(A6 - He) + cos h cos Q^q sin(A6 - Re)) 

(24) 

and 



z,; = - sin/e sin {\q - 



(25) 



where /g, fZe ^-nd Ag are respectively Titan's mean inch- 
nation, argument of the node and mean longitude in the 
inertial frame of the ephemerides. (The subscript 6 refers 
to the fact that Titan is Saturn's sixth satellite.) 

To obtain the coordinates x, y, and z of Saturn is the 
reference frame bound to Titan (/i, /2, /a), 5 rotations are 
to be performed: 



with 



R3{-l)Rii-J)R3{-g)Rii-K)R3i~h) 



(cos (j) — sin (j) ' 
sin (j) cos (f) 
1, 




and 



RM 



10 

COS (j) — sin ( 
sindi cos( 



(26) 



(27) 



(28) 



Table [T] gives the values of the physical and dynami- 
cal parameter that we use, and Table [5] gathers the com- 
puted values of the corresponding parameters used in the 
Hamiltonian (PT|) . 



Table 2. Values used in the Hamiltonian (PT|) that have 
been computed from the physical and orbital parameters 
given Table [TJ 



Parameter 


Numerical value 


do 


1222345.284 km 


n* 


143.8339397847rad.y-i 


71 


1.016129 X lO""* 


72 


7.248387 x 10"^ 


5i 


-1.522286 X 10-* 




-1.085897 X 10"* 


Ss 


9.247193 X 10"^ 



3.1. Equilibrium 

We consider here that the system is exactly at the Cassini 
state. This implies that: 



The ascending node of the frame (ni, ^2, na) (associ- 
ated to the angular momentum) in the inertial frame 
has the same precession rate as the ascending node of 
Saturn in the same inertial frame, i.e. r + f2|^ = 0, 
is the argument of the ascending node of Saturn. 
There is no wobble, so the angular momentum is colin- 
ear with Titan's axis of highest inertia f^. This implies 
J = 0, so = and r]q = 0. 



We have 



and 



6: 



(29) 



(30) 



so it is convenient to introduce this new set of canonical 
variables: 



P ■ 



P 
R 



where a represents the angle between the axis of least iner- 
tia of Titan /i and the direction Saturn-Titan, and p is the 
difference between the two ascending nodes. At the exact 
equilibrium, these two angles should be zero. 

This way, and also assuming d ~ do (i.e., neglecting 
Titan's orbital eccentricity) , the Hamiltonian (PTjl becomes 



n 



nP + ^R 



(31) 



Ss)[ai sin^ K + a2 sin K cos K cos p 
+a3 cos 2p(l — cos 2K)] 
+nS2{l + 5s)[bi{l + cos if) ^ cos 2(7 
+b2 sin if (1 -I- cos if) cos(2cr + p) 
+b3 sin^ K cos(2(T + 2p) 
+b4 sin if (1 — cos if) cos(2(t + 3p) 
+65(1 - cos if )^ cos(2ct + 4p)], 



with the mean longitude disappearing after averaging, ex- 
cept of course in the p variable. The term —nP + flR has 
to be added because the canonical transformation we use 
is time-dependent. The Hamiltonian (j3ip has been com- 
puted with Maple software, and the analytical expressions 
of the coefficients Ui and bi are the same as in Henrard 
and Schwanen (|2004p . This means that, assuming that the 
ascending node of the orbit of the considered body circu- 
lates or not does not change the expressions of Ui and bt, 
the formulae given in Henrard and Schwanen p004p can be 
applied to bodies whose node does not circulate, e.g. J-4 
Callisto, S-6 Titan or S-8 lapetus. The analytical expres- 
sions of these coefficients are recalled in Appl^ while Table 
|3] gives their numerical values in our context. 



1^ = 0, and ^ = = 0- These two last equations give 



At the exact equilibrium, we have tr = 0, p = 0, ^ 



dp 



an 



The axis of least inertia, /i, 
mass of Saturn, so we have p - 



points to the center of Ei{P, K) ^ n P - 1 + (l + 5^) A 



A 







longitude of Saturn in the frame (/i, /2, /a) 



A-j^ as the mean 



cos if — 1 



P sin if 



= (32) 



and 
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Table 1. Physical and dynamical parameters. 



Parameters 



n 
e 

7 ■- 
R 



M 
R 

J2 

C22 

C 



Values 



References 



143.9240478491399rad.y~ 
0.0289 

5.6024 X 10"^ 
58232 km 

1.6298 X 10-2 
2.36638 X 10-*Mi 
2575 km 
3.77747586645 x lO^^fcm^j/" 
(3.15 ±0.32) X 10-5 
(1.1235 ±0.0061) X 10^5 
0.31 



TASS1.6T995] 
TASS1.6 1995] 
TASS1.6 1995] 
lAU 2000 \1MI[ 
Pioneer & Voyager [1959] 
Pioneer & Vovager ll989l 
lAU 2000 12002] 

Pioneer, Voyager + lERS 2003 
Cassini '2006] 
Cassini i200^ 
(...) 



Note 1. The mean values of Titan's mean motion n, eccentricity e and inclination 7 come from TASS1.6 theory (Vienne & Duriez 
I1995|) . the radii come from the lAU 2000 recommendations (Seidelmann et al. I2002p . Titan's mass M and Saturn's J2 come from 
the Pioneer and Voyager space missions (Campbell & Anderson 1989). These two values are those used in TASS1.6 theory, we 
choose to keep them in order to remain coherent. The mass of Saturn has been derived from the fly-bys of the Pioneer and Voyager 
space missions, but the published value is given in solar masses. That is why we also indicate lERS 2003 as a reference, which 
gives us the solar mass. Titan's J2 and C22 come from the fly-by Til of the Cassini space mission (Tortora et al. I2006p . but 
unfortunately no value for j^j^ is available yet. We can only hypothesize that it should be included between 0.3 and 0.4, as the 
case for the Galilean satellites of Jupiter. 



Table 3. Numerical values of a, and 6, 



Parameter 


Numerical Value 


ai 


-4.9990584229813 x 10 


-1 


02 


1.12039208002146 x 10" 


2 


as 


1.56929503117011 x 10" 


5 


bi 


2.49984306803404 x 10" 


1 


62 


5.60213623923651 x 10" 


3 


63 


4.70788509351034 x 10" 


5 


64 


1.75839129195533 x 10" 


7 


65 


2.46284149427823 x 10" 


10 



nA 



PsinK 



(33) 



with 



A = Si [ai sin 2K + 02 cos 2K + 2a^ sin 2K] 
+(52 [ - 2hi sinX(l ± cos if) ± 62 (cos if ± cos2ii:) 
±^3 sin 2K ± 64 (cos K - cos 2K) 
+265 si\iK{l- conK)]. 



(34) 



Since Titan's ascending node oscillates around a fixed 
value, we have f^e = 0. A numerical resolution of (|32]) and 
gives 



K* = 1.1204858615 x IQ-^rad 
= 2311.168arcsec= 38'31.168" 

P* = 1; 

hence, 

R* = 6.2773771522 x 10"^ 
the asterisk meaning "at the equilibrium". 



(35) 
(36) 

(37) 



In the orbital model we use, Jg is constant at 
0.011204858615 rad, which is exactly the value of K* we 
get. Such an accuracy of 11 digits is given to indicate the 
numerical equality of the two values, but it has no real phys- 
ical meaning, so Titan's mean obliquity (measured with re- 
spect to its orbital inclination) should be nearly zero. This is 
confirmed by this formula, given by Henrard and Schwanen 
(|2004]) : 



K* 



Si 



(38) 



which becomes K* « / when the mean value of the preces- 
sion rate of the line of nodes is zero. 



3.2. The fundamental frequencies of the free librations 

Since the equilibrium has been found, the Hamiltonian is 
centered in order to study the behavior of the system near 
the equilibrium. We introduce a new set of canonical vari- 
ables: 



ip^ P 



77, = P - F* 
r]p = R- R* 

Vq 



As a translation, this transformation is canonical. In these 
variables, the main part of the Hamiltonian of the problem 
is quadratic. Its quadratic part is named J\f and we have 



n{l + Ss) 



(39) 



P-qqVq ■ 



The analytical expressions of the coefficients fj,xx and 
jxx are recalled in ApplBj and their numerical values are 
gathered in Table U] The reader should be aware that these 
coefficients are similar to those in Henrard and Schwanen 
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Table 4. Numerical values of the coefficients ^^x and Table 5. Numerical values of the coefficients of the 



Parameter 


Numerical Value 






7<rCT 


2.1717941364 x 10"* 


Parameter 


Numerical Value 


lap 


1.3632529077 x 10"* 


Ci 


2.1717941364 x 10"* 


Ipp 


1.6372888272 x 10"* 


C2 


1.6372032547 x 10"* 


111 


3.4788181236 x 10"* 


C3 


3.4788181236 x 10"" 




5.0000000409 x 10"^ 


i>i 


0.5000000000 


flap 


-6.5206613577 x 10"^ 


V'2 


1.0387557096 


Upp 


1.0387557096 


V'3 


1.4564940392 x 10"^ 


flqq 


1.4564940392 x 10"^ 


a 


-6.1404734778 x 10"'' 






p 


6.2770815833 x 10"^ 



(2004), but different from the ones used by Henrard for lo 
(2005a) and Europa (2005b) where other variables are used. 

We now introduce the following new set of canonical 
variables: 



£_p = axi + (1 — aP)x2 



77cr = (1 - af3)yi - ay2 
Vp = Pvi + 2/2 



with a and (3 conveniently chosen so as to untangle the 
variables ^ and r]. ft can be easily checked that this trans- 
formation is canonical, because it preserves the differential 
form; i.e. 



d£,a-ri'y + d(,p.rip + d^q.r]q = dxi.yi + dx2.y2 + dx^.y^. (40) 

With these new variables, the Hamiltonian ([M]) can be writ- 
ten as 



with 



Cixf + (2x1 + Csxl + ipiyl + ip2yl + i^ivl (41) 



Ci = 7<T<T + 27<,pa + -ippO^ (42) 



C2 = 7<.<./3' - 2/3(1 - a/3)7-p + 7pp(1 - «/3)' (43) 



i^i = i^aail - aP'f + 2/3(1 - a/3)/i,p + /^Vpp (44) 



V'2 = a^/^crcr - 2a/lcp -I- flpp 



1p3 



(45) 



(46) 



(47) 



The numerical values of these coefficients are gathered 
Table El 

We can now introduce the last following set of polar 
canonical coordinates: 



xi = \/2UU* sin? 



X2 = V2VV*smv 




with 



U* 



V* 



w* 



(48) 
(49) 
(50) 



The purpose of this last canonical transformation is to 
show the free librations around the exact Cassini state. The 
arguments of these free librations are u, v, and w, and the 
amplitudes associated are proportional to VU, VV, and 
\/W respectively. We can easily check that this transforma- 
tion is canonical because we have du.U + dv.V + dw.W = 
dxi.yi + dx2-y2 + dx^.y^. We can now write 



N = LOuU + LOyV + UJyjW 



with 



W„ = 2^6 vV^(l + <5s) 
= 2n6\/0^(l -I- 5s) 



Ss) 



(51) 

(52) 
(53) 
(54) 



uj^ 2n6VV^(l ■ 

The numerical results are gathered in Table [H so the 
periods associated to the 3 free librations around the equi- 
librium state are respectively 2.09, 167.37, and 306.62 years. 
Table [7] gives an application of the formulae given in this 
paper to the Galilean satellites of Jupiter lo and Europa, 
and we make a comparison with the analytical results of 
Henrard (|2005al and 2005c) and the numerical results of 
Rambaux and Henrard (2005) obtained with the SONYR 
model (Rambaux & Bois I2004j) . which is a relativistic N- 
body model. The small differences between our results and 
Hcnrard's analytical results come from Henrard neglecting 
03 J ^3 J ^4, and 65 for lo, and 64 and 65 for Europa. 

Table 6. The free librations around the equilibrium state. 



Proper Modes 


uj {rad.y ^) 


T (period in years) 


u 


2.9998383244 


2.0945079794 


V 


3.7541492157 x 10"^ 


167.36642435 


w 


2.0491499350 x 10"^ 


306.62399075 
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Table 7. Comparison of the periods of the free hbrations 
of lo and Europa given by different models. 



Proper Modes 


Henrard 


SON YH 


this paper 


lo 








u 


13.25 days 


13.18 days 


13.31 days 


V 


159.39 days 


157.66 days 


160.20 days 


w 


229.85 days 




228.53 days 


Europa 








u 


52.70 days 


55.39 days 


52.98 days 


V 


3.60 years 


4.01 years 


3.65 years 


w 


4.84 years 




4.86 years 



Note 2. The results labelled "Henrard" come from (Henrard 
I2005a[l for lo and (Henrard 2005c ) for Europa, while the resuHs 
labelled "SONYR" come from (Rambaux & Henrard l^nPS]) . 



Table 8. Initial conditions chosen for the numerical inte- 
gration, at t=-4500 years. 



Variable 


Expression 


Po 


Aso — 71" 




ro 






Co 


10"* 




■no 


10-* 




Po 




(-10-* 


Ro 


1.0001 X Po(l - cos 


Ko) 



Note 3. These conditions have been arbitrarily chosen near the 
Cassini state, with Aso, fieoi fieo respectively the values of 
Titan's mean longitude, argument of the ascending node, and its 
instantaneous angular rate, at t=-4500 years, given by TASS1.6. 
Ko is the initial value of the obliquity K on the same date. 



4. Numerical study 

To check the reliability of our previous results and to go 
further in the study of Titan's forced rotation, we used the 
numerical tool. This allowed us first to obtain a solution for 
the rotation of Titan and then to describe it by frequency 
analysis and to split the free from the forced solutions. 

4.1. Numerical integration 

We integrated the 6 equations coming from the Hamiltonian 
(PT|) over 9000 years, i.e. between -4500 and 4500 years, the 
time origin being J1980. In these equations, x and y come 
from TASS1.6 ephemerides. We recall that x, y, and z are 
the coordinates of the barycenter of mass of Saturn in the 
frame (/i, /2, /a) rigidly linked to Titan (see Sect. 2.3). 

We used the Adams-Bashforth-Moulton 10*''-order 
predictor-corrector integrator, with a constant timestep 
/i = 1.6 X 10^* year, i.e. 5.844 x 10^^ day. We considered 
that the shortest significant fundamental period of the sys- 
tem is given by SAe, i.e. ~ 5.315 days w 90 x /i. 

Table |8] gathers the initial conditions we used. These 
conditions were arbitrarily chosen near the Cassini state. It 
implies that, by choosing these initial conditions, we sup- 
posed that Titan is at the Cassini state. In these initial 
conditions, the initial value of K Kq is defined as 



This equation is very similar to ([55)1 . We used it to be sure 
that the system is near the equilibrium. We did not want to 
start at the exact equilibrium but very close in order to be 
able to detect the 3 free librations that we studied in the 
previous section. However, we should keep in mind that 
the frequencies computed in Table [6] are in fact limits of 
the frequencies of the free librations when their amplitudes 
tend to zero. Thus, too high amplitudes of free librations 
would alter the frequencies too much. In that way, their 
comparison with the expected fundamental frequencies of 
the free librations might be difficult, so their identification 
as these fundamental frequencies could become doubtful. 

Figure [2] gives plots of some significant data resulting 
from the numerical integration. Figure [5^ shows the be- 
havior of the variables P = ^ (modulus of the plotted 
value) and a, and Figure [2]d shows the behavior of p, i.e. 
the difference between the two nodes. We can see that this 
angle is oscillating around 0, as predicted by the theory. We 
can also visually detect a period of about 700 years, and will 
see later that it is a forced component due to the behavior 
of Titan's orbital ascending node. Figure O; shows the be- 
havior of the wobble J, we obtained it from the variables 
and rjq. Finally, Figure [2li shows the "obliquity" K . In this 
last panel, we can see the same 700-year-periodic contri- 
bution detected in p. Unfortunately, looking at these plots 
does not give information on the free and the forced com- 
ponents of the solutions. That is why we used the frequency 
analysis technique. 

4.2. Analysis of the solutions 

We use the frequency analysis to describe the solutions 
given by the numerical integration, i.e. to give a quasi- 
periodic representation of these solutions. Such a technique 
has already been used often to describe the orbital motion 
of planets (see Laskar ll988p or natural satellites (cf. for in- 
stance Vienne & Duriez 119951 or very recently Lainey et al. 

mm . 

One of the main difficulties with this kind of problem 
is that we have two timescales for the periods of the terms 
that appear in the synthetic representations: Titan's orbital 
period is roughly 16 days, while the period of its pericenter 
is about 700 years. In order to correctly detect the long- 
period terms, the total time-interval used to analyze the 
solution should be about as long as the longest period ex- 
pected, here « 3200 years, and the timestep shorter than 
half the shortest period expected (about 5 days). Thus, data 
over 3200 years should be represented with a timestep of 
2.5 days, but this would require about 500000 points. This 
would take a very long computation time, but fortunately 
some alternative techniques exist to solve this problem. 

The most common technique is the use of a digital filter 
that splits the short-period terms from the long-period ones 
(see Carpino et al. I1987p . However, this technique might 
alter the signal. Another technique has been used here, 
which consists in using two samples of data with very close 
timesteps, as explained in Laskar l|2004p . More precisely, 
for each variable, we extracted two samples of 65536 data 
from the results of the numerical integration, one point ev- 
ery 848 for the first sample and one point every 864 for the 
other one. As a result, the first sample represents the so- 
lutions for 8891.7888 years with a timestep hi = 49.55712 
days, and the second represents the solutions for 9059.5584 
years with a timestep /12 — 50.49216 days. 
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time (kvr) 
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time (l<vr) 
(d) 

Figure 2. Numerical simulation of Titan's obliquity over 9000 years, the time origin being J1980= 2444240 JD. Here 
the behavior of the variables P, ti, p, J and K is being displayed. Explanations are in the text. 



These two timesteps are far too large to detect contri- 
butions with a period of about 16 days. In fact, the short- 
period terms are detected, but with a wrong frequency. 
When a frequency v is too high, it is detected as 



vi = — 
hi 

in analyzing the first sample, and as 



V2 



h2 



(56) 



(57) 



in analyzing the second one, where ki and ^2 are (a priori 
unknown) integers. We have 



{V2 - Vl)h2 ^ k2~ ki — 
hi 



(58) 



and 



ki. 



(59) 



V2h2 - vihi = v{h2 - hi) + k2 

If we now define [x] as the closest integer to the real x (i.e. 
\[x\ — x\ < i), we have 



[v2h2 - i^ihi] = k2- ki 



(60) 



and finally 



fci = 



-((j/2 - l^l)h2 - [j^2/l2 - l^lhl]), 



(61) 



hi — /i2 

where Eq. (l60l) requires that hi and /12 are close enough, i.e. 
|z/(/i2 — hi)\ < 5- In our case, the highest frequency that we 



can detect with this method is 



0.54c? , so we 



2(/l2-/ll) 

can detect every term with a period longer than 11.75 days, 
while analyzing only one sample would give periods longer 
than about 100 days (i.e. 2 timesteps). Such accuracy is 
enough to detect Titan's orbital period. 

Table [9] is an example of the decomposition of a variable 
(here P = with G the norm of the angular momentum) 
with the two timesteps. The algorithm used for determin- 
ing each frequency is taken from (Laskar et al. I1992P and 
has been iteratively applied to refine each frequency, as de- 
scribed in fChampenois I1995|) . In this table, term 1 is a 
constant part, while the second one is clearly the free libra- 
tion associated to the proper mode u. The slight difference 
between the obtained and the expected periods should be 
partly due to the associated amplitude not being null, and 
partly to the approximations used in our analytical model 
(i.e. no eccentricity and a constant inclination). However, 
we can see that the two determinations give very differ- 
ent results for term 3, so we can infer that it is in fact a 
short-period term. 
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Table 9. Decomposition of the solution for P with the two timesteps, i.e. hi = 49.55712 days (left) and /i2 = 50.49216 
days (right). 



N° 


Amp. 


Phase (°) 


T(y) 


Amp. 


Phase (°) 


T(y) 


1 


1.000000002 


-4.57 X 10"** 


3.50 X 10" 


1.000000002 


8.95 X 10""* 


-1.82 X 10" 


2 


0.000099514 


0.76 


2.09773 


0.000099514 


0.77 


2.09773 


3 


0.000025104 


144.00 


1.25952 


0.000025104 


144.00 


0.83097 



Note 4- The origin of phases is here the origin of the frequency analysis, i.e. 4499.99344 years before J1980. The series are given 
in cosine. 



Applying (|6ip we find T = 15.6612 days. That is very 
near to Titan's orbital period, so this term should be an in- 
teger combination of Titan's orbital period and other con- 
tribution(s). 

It is now interesting to identify the periodic terms con- 
tained in the solutions associated to the considered vari- 
ables. Table [TU] gives the proper modes that are expected. 
They should appear in the quasiperiodic decompositions of 
the solutions as parts of integer combinations, so integer 
combinations of the frequencies of the proper modes are 
performed to identify each term of the decompositions. We 
do not use the phases because they are uncertain in the 
Titan ephemerides given by Vienne & Duriez |1995j). The 
reason is that the given phases are in fact integer combina- 
tions of the phases coming from the identified proper modes 
and very-long-period arguments due to the solar perturba- 
tion that are assumed to be constant on an ephemerides- 
timescale. 



terms. Except for two of them, all the components have 
been clearly identified. The terms whose periods are written 
in days were determined by comparing the results given by 
the two analysis. These terms are not clearly identified, the 
reason probably being that they require a high accuracy in 
their determination. In (|6ip . two quantities are substracted, 
so a cancellation problem might appear and complicate the 
determination. Moreover, an integer combination between 
a short-period term (like Titan's mean longitude Ag) and 
a long-period term gives a short-period term very close to 
the original short period, so the short-period terms that we 
detected might in fact be sums of several terms with very 
close frequencies, making them very difhcult to split. 

These short-period terms seem to have a period very 
close to Titan's orbital period, except the term 5 in the 
decomposition of a. If we consider that the timesteps hi 
and /i2 are not close enough and that we have in fact [i'{h2 — 
hi)] = 1, (|6T|) becomes 



Table 10. Proper modes of the system. 



Proper 


Frequency 


Period 


Cause 


Mode 


[rad.y^^) 






As 


508.00932017 


4.52 days 


Rhea 


Ae 


143.92404729 


15.95 days 


Titan 


As 


28.92852233 


79.33 days 


lapetus 


4>5 


0.17554922 


35.79 years 


65 


$5 


-0.17546762 


35.81 years 


75 


(t>6 


0.00893386 


703.30 years 


66 


$6 


-0.00893124 


703.51 years 


76 


<Pa 


0.00197469 


3181.86 years 


68 


$8 


-0.00192554 


3263.07 years 


78 


As 


0.21329912 


29.46 years 


Sun 




2.995 


2.09773 years 


Vu 


4>v 


0.0375 


167.4883 years 


W 


4>m 


0.0205 


306.3360 years 


Vw 



Note 5. The modes A5 to A9 (first part of the Table) come from 
Vienne & Duriez (|1995p . while the second part contains the free 
librations around the Cassini state. These terms have been eval- 
uated from the solutions given by our numerical integration. The 
fourth column gives the orbital parameter to which the proper 
mode is linked, being the eccentricity of the satellite i, and 7^ 
the sine of its semiinclination. The subscripts i are 5 for Rhea, 
6 for Titan, and 8 for lapetus. 



The results are summarized in Tables [Tl] to [T6l the ori- 
gin of the phases now being J1980. There, K is given in 
radians, and the other variables have no unit. The terms 
with a period T written in years could have been obtained 
in analyzing only one set of data. In fact, the two sets 
have been analyzed, and the results are the same for these 



'^2 

/Ci = 1 + T— ((i^2 - l^l)h2 - [V2h2 - vihi]), (62) 

hi — ri2 

and we obtain a term whose period is 5.22008 days. This is 
quite close to the period associated to 3Xq. 

The quasiperiodic decomposition of the solutions allows 
us to split the forced solution away from the free one. The 
free solution around the equilibrium can only be known 
with observations that could give initial conditions for the 
numerical integration. However, the forced solution only de- 
pends on the equilibrium and can be obtained in dropping, 
in the solutions given in Tables [TT] to [TBI the terms depend- 
ing on the free libration modes 4>y and 0^,. Thus, we can 
for instance see that angle J is not zero at the equilibrium 
but has a forced motion. This possibility of a forced wobble 
has already been pointed out by Bouquillon et al. (|2003p in 
a general study of the rotation of the synchronous bodies 
(i.e. that are in a 1 : 1 spin-orbit resonance). 

5. Discussion 

5.1. Comparison between the analytical and the numerical 
results 

Table [T7] gives a comparison between our analytical and 
numerical results. We recall that, in the analytical model, 
the orbit of Titan is circular with a constant inclination, 
whereas the orbital eccentricity of Titan (i.e. 0.0289) is 
taken into account in the numerical model, along with the 
variation in its inclination. We can see very good matching 
for the periods of the free librations around the equilibrium. 
In contrast, we can see a significant difference in the equi- 
librium obliquity K*. The line e refers to the equilibrium 
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Table 11. Quasiperiodic decomposition of the variable P. 



N° 


Amp. 


Phase (°) 


T(y) 


Ident. 


Cause 


1 


1.000000002 


4.89 X 10"'" 


-1.82 X 10"'' 


constant 




2 


0.000099514 


63.00 


2.09773 


(t>u 


Vu 


3 


0.000025104 


35.44 


15.6612 days 


unknown 





Note 6. The series are in cosine. The fourth column gives the orbital parameters associated to each identified term. 



Table 12. Quasiperiodic decomposition of the variable R. The series are in cosine. 



N° 


Amp. xlO® 


Phase (°) 


T(y) 


Ident. 


Cause 


1 


9.17912502 


2.07 X 10"" 


6.12 X 10'" 


constant 




2 


8.22242693 


-170.93 


703.50790 


-$6 


76 


3 


2.28952572 


174.72 


167.48834 


0t) 


Vv 


4 


1.49857669 


-143.60 


219.82166 


$6 + (t>v 


W 76 


5 


0.30469561 


-107.54 


3252.81 


$8 


78 


6 


0.22920477 


-61.85 


899.49195 


$8 - 


7678 


7 


0.05732146 


-80.14 


176.5406 


$8 + (Pv 


V^78 


8 


0.01907471 


163.85 


14.72857 


2A9 


Sun 



Table 13. Quasiperiodic decomposition of the complex variable riq + \/—l^q. 



N° 


Amp. xlO" 


Phase (°) 


T(y) 


Ident. 


Cause 


1 


9.12391728 


-51.69 


306.33602 




y/W 


2 


6.01688587 


51.69 


-306.33605 


—4>w 


Vw 


3 


5.73033451 


158.48 


351.70284 


06 - 


6676 


4 


3.83212940 


-158.48 


-351.70284 


<E>6 - 06 


6676 


5 


0.63642954 


-35.86 


135.27368 


0„ - $6 




6 


0.38395548 


35.86 


-135.27368 


$6 - (t>v 


W 76 



obliquity with the normal of Titan's orbit as its origin. It is 
computed by substracting the mean inclination of Titan to 
K*. The mean inchnation of Titan is 1.12049 x 10^^ in the 
analytical model and 1.18985 x 10^^ in the numerical one. 
The difference in K* partly comes from the difference in 
the mean inclination of Titan, but probably not only from 
it. 



5.2. Influence of Titan's inclination and eccentricity 

Titan's inclination plays an overwhelming role in its obliq- 
uity, as shown in Moreover, the proper modes $5, ^e, 
and $8 given in Table [10] can be linked in its inclination, 
because they consist of the main (or at least second) part 
of the solutions for C = sin |ea;p(v'— 1^2) for Rhea, Titan, 
and lapetus in TASS1.6, and they appear in the solution 
for Ce (related to Titan's inclination). It is striking, for in- 
stance in reading Table [151 that the term $6 plays a very 
important role in the forced and in the free solution. We 
can even figure the period of 703.51 years just in looking at 
Fig[llD. 

In contrast, the proper modes (^5, (/>6, and 08 do not 
clearly appear, with the exception of 06 in rjq + 
These modes are related to the eccentricities of Rhea, Titan, 
and lapetus, and their values are respectively lO^'^, 0.0289, 
and 0.0294. In fact, 06 might play a more important role 
than suggested by Tables [TT] to [TBI because it could be con- 
fused with — $6 by the algorithm of frequency analysis. The 
reason is that these two terms have very close periods, i.e. 
703.3 and 703.51 years. If we call vi and 1^2 the associated 
frequencies and vq the Fourier fundamental frequencies (i.e. 



the frequency associated to a term whose period is the in- 
terval of study, 9000 years in our cases, it has no link to 
the fundamental frequencies of the system), the algorithm 
of frequency analysis can split I'l from 1^2 only if 

-z/2| > 21/0- (63) 

This implies that the interval of study should be longer than 
4.712 X lO^yr. Such a timcscale is not consistent with the 
ephemerides and so cannot be considered. It does not mean 
that the terms identified as $6 are in fact 06 1 because the 
period found is much closer to 703.51 years than to 703.30. 
It just means that there might be a very small contribution 
due to 06 in the identified term. Moreover, we cannot ex- 
clude a role played by the eccentricities in the values of the 
equilibrium obliquity and of the fundamental frequencies of 
the free librations. 

5.3. Uncertainty on Titan's gravitational field 

Titan's gravitational field is not clearly known. We are con- 
fident in its mass thanks to the Pioneer and Voyager fly- 
bys (see Campbell & Anderson ll989p . but we are uncertain 
about 10% of its J2, and we have no value for the ratio 
jY]^- We can just hypothesize that it is included between 
0.3 and 0.4, as it is the case for the Galilean satellites for 
Jupiter. We arbitrarily chose jf^ = 0.31 and also tried 
with jjj^ = 0.35. 

Tables [TS| and [Hj and Figure [3] summarize the result of 
the study of Titan's rotation with jjj^ = 0.35. Except for 
(c) , the plots do not show any evident difference with Figure 
[21 because the frequencies of the free librations are shifted 
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Table 14. Quasiperiodic decomposition of the variable a. The series are in sine. 



N° 


Amp. xlO^ 


Phase (°) 


T (y) 


Ident. 


Cause 


1 


4.78176461 


63.00 


2.09773 




Vu 


2 


0.02510524 


35.43 


15.6612 days 


unknown 




3 


0.01147635 


-5.61 


167.47831 




Vv 


4 


0.01014094 


8.48 


703.51797 




76 


5 


0.00961527 


-74.09 


55.1128 days 


unknown 




6 


0.00896290 


-56.74 


2.11032 


(t>n + 2$6 




7 


0.00891238 


2.47 


2.08529 


<t>u - 2$6 




8 


0.00744233 


137.13 


15.6602 days 


unknown 




9 


0.00566630 


165.69 


219.75766 


(t)v + "1>6 


W 76 



Table 15. Quasiperiodic decomposition of the variable p. The series are in sine. 



N° 




AlTlp. 


Phase (°] 


T (v) 


Idsnt . 




1 





1 or^onooT 

18U898o7 


175.64 


167.49723 




/t7 
V V 


z 


U 


loDD / ooy 


1 vn on 


/ Uo.0z44D 


rh 

— cpg 


76 


3 





11829380 


— 175.17 


IOC O OTO 4 

135.28724 


J. /Fx 

<Pv — '±'6 


V 1/76 

„2 
76 


4 


U 


09023900 


— iOi.9i 


OKI TKTOn 

3oi. (0 (89 


2'I>6 


5 





07735641 


-166.02 


113.46712 


4>v - 2$6 




6 





05226443 


-152.60 


234.50407 


-3$6 


76 


7 





05058400 


-156.89 


97.70793 


(t)v - 3$6 


/T7 3 

VV^76 

76 


8 





03311443 


-147.68 


85.79329 


(pv - 4$6 


9 





03060799 


-143.30 


175.88361 


-4$6 


10 





02165111 


-138.40 


76.46679 


4>v - 5<E>6 


VV^76 


11 





02035272 


-173.58 


74.84659 


2(j}y - $6 




12 





02027563 


-167.25 


67.64828 


2(j>^ - 2$6 


^7l 


13 





01799820 


-132.60 


140.70722 


-5$6 


7I 


14 





01785634 


-158.19 


61.71406 


(t)v - 3$6 


^76 


15 





01739585 


163.87 


14.72858 


2A9 


Sun 


16 





01498249 


-176.87 


83.76245 


2(j}v 


V 


17 





01472909 


-149.81 


56.73966 


2(t>^ - 4*6 


Vlt 


18 





01417957 


-128.79 


68.97041 


(j)v - 6$6 




19 





01152268 


-140.62 


52.50108 


2cj>^ - 5$6 


20 





01065513 


-121.56 


117.25730 


-6$6 


7I 


21 





00929051 


-148.93 


62.81242 


unknown 




22 





00899501 


154.97 


15.04352 


2A9 + $6 


Sun, 76 


23 





00881939 


-131.88 


48.85511 


2<j>^ - 6*6 


Vll 


24 





00657942 


-122.55 


45.68121 


2<l>y - 7$6 


Vll 


25 





00635299 


-110.18 


100.50834 


-7$6 




26 





00609640 


-109.01 


57.66387 


unknown 




27 





00584381 


-87.12 


109.64279 


(t)v - 2<E>6 - <E>8 


^76 78 


28 





00584403 


-95.53 


129.88963 




7678 


29 





00531164 


-19.68 


29.44635 


A9 


Sun 


30 





00513296 


-81.36 


317.28640 


-2$6 - $8 


76 78 



Table 16. Quasiperiodic decomposition of K. The series are in cosine. 



N° 


Amp. xlO^ (rad) 


Phase (°) 


T(y) 


Ident. 


Cause 


1 


1.25481164 


8.68 X 10""' 


-2.65 X 10" 


constant 




2 


0.68465799 


-170.92 


703.51272 


-$6 


76 


3 


0.17842225 


175.02 


167.49146 






4 


0.10246867 


-161.88 


351.76856 


-2<E>6 


5 


0.07264971 


-15.67 


219.80041 


(i>v + "I>6 


^A^76 



just a little when -^j^ changes. In contrast, the behavior 
of the wobble J (Figure[3t) is very interesting, because this 
angle can be 10 times bigger than in the previous simula- 
tion. Table [T51 indicates that the most important terms in 
the solutions of and 77^, on which J depends, are and 
2$6- In our cases, the periods of these terms are very close, 
so there might be a resonance between them, which could 
explain the amplitude of J. The matching on the frequen- 



cies of the free librations between the analytical and the 
numerical methods is still good, while a shift on Titan's 
mean obliquity still exists. 

6. Conclusion 

This paper offers a first study of Titan's rotation, where 
Titan is seen as a rigid body. We obtain a quasiperiodic de- 
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Table 17. Comparison between our analytical and numerical results. 



Parameter 


Analytical 


Numerical 


Difference 


K* (rad) 


1.1204859 X 10~^ 


1.25481164 X 10"'' 


12% 


e (arcmin) 





2.233 


(...) 


Tu (y) 


2.094508 


2.09773 


0.15% 


n (y) 


167.36642 


167.49723 


0.08% 


(y) 


306.62399 


306.33602 


0.09% 




time (kvr) 
(c) 




time (l<vr) 
(b) 




time (l<vr) 

(d) 



Figures. Numerical simulation of Titan's obliquity over 9000 years, with — 0.35. The displayed variables are the 
same as in Fig (31 



Table 19. Comparison between our analytical and numerical results, with 



c 



= 0.35. 



Parameter 


Analytical 


Numerical 


Difference 


K* (rad) 


1.1204859 X 10"^ 


1.272996 X 10~^ 


13.6% 


e (arcmin) 





2.858 


(...) 


Tu (y) 


2.225839 


2.22896 


0.14% 


Tv (y) 


188.987571 


189.10854 


0.06% 


Tu, (y) 


346.236493 


348.49661 


0.65% 



composition of the forced solution, which can be split from 
the free solution in which Titan's obliquity plays an over- 
whelming role. Moreover, we find good matching between 
the frequencies of the free librations around the equilibrium, 
analytically and numerically evaluated, despite a model of 
circular orbit in the analytical study. However, we find a 
slight difference in the equilibrium obliquity. Finally, we 



cannot exclude a resonance between the proper mode $6 
and Titan's wobble. 

The next fly-bys of Cassini spacecraft should give us 
more information on Titan's gravitational field, so we 
should be able to make a more accurate study on its ro- 
tation, that could include direct perturbations on the other 
Saturnian satellites. These perturbations are supposed to 
be small (see for instance Henrard l2004p and should be neg- 
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Table 18. The free librations around the equiUbrium state, 
with = 0.35 



Proper Modes 


w {rad.y ^) 


T (period in years) 


u 


2.822839 


2.225839 


V 


3.324655 x 10"^ 


188.987571 


w 


1.814709 X 10"^ 


346.236493 



ligeable compared to the uncertainties we have on Titan's 
gravitational parameters. After that, the next step is to 
consider Titan as a multilayer non-rigid body and to study 
the consequences of its internal dissipation on the rotation. 
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Appendix A: The coefficients a^, 6j 



ai = 



sin'' / 1 + cos^ / 



a2 = ^ = 27^1^(1 - 27^) (A.2) 
a3 = ^ = ^(l-7^) (A.3) 



b2 



_ 1 2 cos/ + cos^ / _ 1 - 272 -h 7^ 



16 

2 sin 7-1- sin 27 



: 7(1-72)71^ 



63 = ^sin2 7=^7^(l-y) 



64 = 



2 sin I — sin 27 



1 — 2 cos I + cos^ I 7^ 
= 16 = T 



(A.4) 

(A.5) 
(A.6) 
(A.7) 
(A.8) 



Appendix B: The coefficients ^xx and fixx- 

7^^ = -2S2{bi{l + cosK*)^ + b2smK*{l+cosK*) 
+63 sin^K* + bi sin7r*(l - cosK*) + 65(1 - cosK*)^) 



(B.l) 
(B.2) 



-y^p = -52(62 smK*{l -\- cosK*) -f 63 sin^ K* 
-1-364 sin ii'*(l - cosTi"*) + 465(1 - cosK*f) 

Ipp = - (^Si sin 2K* + 4a3 sin^ K*^ 

+52^^ sin7C*(l -h cosTT*) -\- 263 sin^ K*-\- (B.3) 

^64 sin if* (1 - cosK*) + 865(1 - cosif*)^)^ 

1 71+72 



799 = 



-{Sl+S2)[ 



/cos{K*-I) 7 



2 1 - 71 - 72 
5 



+ -cos(2(if*-7)) + -) 

(B.4) 



A*'-'- = \ + +2a3)(l -cos7sr*)(3cos7r* - 1) 

^1 (1 tl^K. f + 4 co.^ - 5 cosif - - 2)) 



(61 (3 cos if* + 1) (cos if* - 1) + -62 



3^ sin if* cos 27^* 
2 ^ l + cos7f* 



-1-63(1 - cos iC*) (3 cos if* - 1) 

^64 1-COS^* ^* ^ 8 ^* ^ g ^Qg2 ^ 3^ ^ 

2 1-1- cos if* / 

(B.5) 

Ma^ = ;^((«i+2a3)(l-2cos7f*)+ 
02 1 -F4cosif* - 2cos2 if* - 4cos3 7f 



2 

P 



+ 



(261 cos if* -62 



sin if* (1 + cos if*) 

4cos2 7f*-cos7f*-2 



2 sin if* 

, T^* .X 64 cos if * cos if * - 1 , ^ , . 

-63 2 cos 7f* - 1 + ^—777 — — 4cosif* + 5 

2 sm if* cos if* + 1 

sin^Tf* \ 
-265: ^ 



1 -|- cosTf * 



(B.6) 
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= -prA^r + y -3^^ + 2«3 

62 /, b2 2cos^K* -3cosK* - 



bil-3cosK* +2cos^ K* 



+ 65) 



^ 2 sin^ K* 

(B.7) 

2 1 - 71 + 72 



